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In this paper, we extend a model of host-parasite co-evolution to incorporate the semi-conservative 
nature of DNA replication for both the host and the parasite. We find that the optimal mutation 
rate for the semi-conservative and conservative hosts converge for realistic genome lengths, thus 
maintaining the admirable agreement between theory and experiment found previously for the con- 
servative model and justifying the conservative approximation in some cases. We demonstrate that, 
while the optimal mutation rate for a conservative and semi-conservative parasite interacting with 
a given immune system is similar to that of a conservative parasite, the properties away from this 
optimum differ significantly. We suspect that this difference, coupled with the requirement that a 
parasite optimize survival in a range of viable hosts, may help explain why semi-conservative viruses 
are known to have significantly lower mutation rates than their conservative counterparts. 

PACS numbers: 87.14Gg, 87.23Kg, 87.10+c 



I. INTRODUCTION 

Introduced over 30 years ago, the quasispecies model 
of evolution^ has provided an invaluable tool for the 
study of complex evolutionary behaviors. In the model, 
a fitness landscape is introduced, which accounts, often in 
a highly approximate manner, for the complex interplay 
between genotype, phenotype and environment by as- 
signing a relative fitness for each genomic sequence (and 
thus associating phenotype with genotype, an approxi- 
mation that must be treated with care). Through the 
consideration of numerous individually mutating copies 
of a genome, evolutionary systems can be studied analyt- 
ically and numerically on these fitness landscapes, which 
has provided enormous insight into the process of evolu- 
tion and the nature of mutation rates in real biological 
systems. In particular, it was found that a phase tran- 
sition (known as the "error catastrophe") occurs as the 
mutation rate increases, and a marked crossover can be 
observed from the existence of a quasispecies (wherein 
most individuals in the population contain genomes close 
to a fitness peak) to a near-random walk in genome space 
with no discernible quasispecies present. 

The vast majority of the literature on the quasispecies 
model involve studies of asymptotic behavior on numer- 
ous stationary landscapesiiSi&. This corresponds to a 
situation where static environmental conditions are con- 
sidered to be the dominant evolutionary pressure on a 
species. However, this picture fails to describe the cor- 
nucopia of evolutionary pressures in nature. Many or- 
ganisms, parasites, survive through the detrimental use 
of host biochemical processes. The parasite requires the 
host to live. The host survives better if it can avoid or de- 
stroy the parasite, providing an intriguing scenario: the 
host must evolve to defeat the parasite and the parasite 
must evolve to evade the host's defenses. This creates a 
non-linear feedback cycle as both species scour a time- 
dependent fitness landscape that changes as the other 
species mutates. 



Parasites are ubiquitous in nature, ranging from the 
microscopic (e.g. viruses, bacteria, protozoa) to fungi, 
helminths and arthropods. The interaction between par- 
asites and hosts are very complex, with parasites exhibit- 
ing multistage life cycles, inert phases, and the use of 
multiple intermediate hosts, while hosts employ a wide 
variety of behavioral and immune defenses. This ongoing 
struggle has been well documented in mammals, birds, 
fish, bacteria and other organisms. 

Recent work on time-dependent quasispecies 
landscapes^ has allowed for the study of a simple 
model of co-evolution by Kamp and Bornholdta— , 
discussed in detail in Section III. They derived a 
parameter-independent expression for the optimal muta- 
tion rate for a host genome, which compared admirably 
with experimental results on B-cell mutation rates^. An 
expression was also derived for optimal viral mutation 
ratesifi which, although dependent on the parameters 
of the model, explained numerous phenomena including 
the constancy of mutation rates within a viral class. 
However, this model considers the interaction only 
between a conservatively replicating parasite and host. 

In its conservative formulation, the quasispecies model 
considers single stranded genomes that produce multiple 
copies of itself, each possessing a set of point mutations, 
while the original genome is conserved. While this model 
is obviously applicable to numerous RNA-based viruses, 
the vast majority of organisms, including many viruses 
and other parasites, store genetic information in dou- 
ble stranded DNA. DNA replicates semi-conservatively 
through a series of steps discussed in Section II. In a re- 
cent work, Tannenbaum et. al. 11 reformulated the qua- 
sispecies model to accurately represent semi-conservative 
systems, which were found to display fundamentally dif- 
ferent behavior than conservative systems with respect to 
the error catastrophe in the infinite time limit on a static 
landscape. Thus, to properly model the co-evolution of 
a parasite and its host, the host system must replicate 
semi-conservatively, while the parasite can be modeled as 
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either conservative, as in the case of many riboviruses, or 
semi-conservative, as by many lysogenic double stranded 
DNA viruses or higher parasites. Retroviruses, such as 
HIV, likely display characteristics of both modes of repli- 
cation, as do immune systems that undergo somatic hy- 
permutation. 

In this paper, we extend Kamp and Bornholdt's model 
of co-evolution to the case of a semi-conservative host in- 
teracting with either a conservative or semi-conservative 
parasite. We consider the optimal behavior for both 
the host and parasite, and demonstrate the similari- 
ties and differences between the conservative and semi- 
conservative models. 

The paper is organized as follows: in Section II we 
present the quasispecies model and its extension to semi- 
conservative replication. In Section III we discuss the 
model of host-parasite co-evolution for both conservative 
and semi-conservative organisms. Section IV presents 
the results and discussion and Section V presents our 
conclusions. 



II. THE QUASISPECIES MODEL 



These equations can be greatly simplified in the case of a 
single fitness peak landscape, where a master sequence, 
<po, has a fitness much greater than all other sequences. 
The rest of the genomes are assumed to be equally fit, 
which can be described by the growth rates 



<J > 77 



(3) 



The sequences can then be grouped into Hamming classes 
based on their distance from the master sequence by 
defining 

wi = x <t> ( 4 ) 

<f>e{4>\HD(4>,4> )=i} 



and 



A{l)=A{<t>) <P^{cp\HD{<P,fa) = l}. (5) 



This reduces the problem from S n dimensions to n + 1 
dimensions. If mutations that lead from higher to lower 
Hamming distances are ignored (an approximation that 
becomes exact as n — > 00), 



In this section, we present some necessary background 
on the conservative and semi-conservative quasispecies 
models for the purpose of a self-contained discussion. 
Greater detail may be found in the original papers. 



A. Conservative Replication 

The quasispecies model studies the evolution of a pop- 
ulation of organisms, each with a genome <fi = S1S2 • • • s n , 
where each represents a "letter" chosen from an al- 
phabet of size S. Often, S is chosen to be two to model 
the pyrimidine and purine groups or four to model the 
nucleotides. Assuming first-order growth kinetics and as- 
sociating phenotype with genotype (i.e. that the growth 
rate of an individual is directly determined by <p) , it can 
be shown that 



(1) 



where denotes the fraction of the population with 
genome <p, A{<j>) represents the fitness, or growth rate, 
of sequence (j), W{(j), </>') is the likelihood of creating se- 
quence <f> from </>' by mutations, and f(t) — A((j))x^, is 
the average fitness of the population, holding the popu- 
lation size constant and introducing competition. If only 
point mutations are allowed and a genome-independent 
mutation probability e is assumed, then W((j>, 4>') can be 
written in terms of the genome length n and the number 
of bases at which <p an d <p' differ, the Hamming distance 
HD(4>,<f/), as 

W(<l>, 4>') = (_f_)*«W)(l - e) n-HD(^')_ (2) 

s 1 



dwi 



y- l (i-ey 



l'=0 



(n-iy. 



-(j-0 



-f(t)wi, 



w v 



(6) 



where f(t) = J2i M l ) w i = aw o + vi 1 - w o) = (o - - 
ri)wQ + r\. Defining y, = iOj exp(J f(s)ds) removes the 
non-linearity in these equations and the linear set of dif- 
ferential equations can be solved for any Hamming class. 
The solution for the master sequence is 



y (t) = y (0)e"^ t 
and, for the first Hamming class, 



(7) 



*w=tt>(°H - {s -i)(A )g J- (8) 

where q = 1 — e, a definition we shall use throughout the 
paper. 



B. Semi-conservative Replication 

In order to properly model a semi-conservative sys- 
tem, a double stranded molecule generated from an al- 
phabet of size S must be considered, where each "let- 
ter" i uniquely pairs with (i + 5/2) mod S. DNA re- 
quires S = 4, where the letters can be assigned as 
A = 1,G = 2,T = 3,C = 4. A single DNA molecule 
of length n consists of a strand <j> = sis 2 ■ ■ ■ s n and a 
complementary strand = ■ ■ ■ s n where s 4 denotes 
the complement of Sj. Hence, each DNA molecule may 
be represented by the pair {<p, cf)} = {</>, </>}. 
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Strand Separation 



produce the pair {<p, </)}. To make these equations more 
useful, we can define ^4(0) = A((f>, </>) and x^ = ^X{tj>,<p} 



if 4> 7^ and x^ 
tion, we obtain 

dx<t, _ 
~df ~ 



if, 



After some manipula- 
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Lesion Repair 



FIG. 1: A schematic model of DNA replication. Adapted 
from Tannenbaum et. alw~. 



When a semi-conservative molecule replicates, it un- 
dergoes a three step process shown schematically in Fig. 
1. First, each genome {(/>, <fi} unzips to form two single 
stranded genomes, 4> and 4>. Each strand is then copied 
to produce two new pairs, {(f>, <fi'} and {</>, (/>'}, where the 
primes denote the fact that the two fresh strands may 
contain replication errors. At this point, proofreading 
mechanisms can distinguish between the new and old 
strands and may fix all or some of the replication er- 
rors, which can be spotted by the fact that ^ s^. 
All of these repair mechanisms are included in the base- 
independent error probability e. In the last step, the 
new and old strands become indistinguishable. Various 
maintenance enzymes repair the remaining mismatches, 
but cannot determine which of the strands <p and 4> is the 
newly replicated strand. Hence, the repair is made in the 
new strand with 50% probability and in the old strand 
with 50% probability. The final result is that the origi- 
nal strand {cf>, <fi} is replicated to create two new strands, 

{<t>",f} and {<t>";,f')XL. 

The quasispecics equations for this system can be writ- 
ten asii 



dx 



{<f>;4>} 



dt 



-- ]T A({<p',^})x {<f ,, >f} (p( ( p',{ ( p,0) + 

p(0', {<M})) - (A({<m» + f(t))x {M , (9) 



where /(f) = J2^, A ({^^)) x {<l>,fl and p{4>\{4>,<^}) rep- 
resents the probability that the unzipped strand 4>' will 



-(A(<t>) + f(t))x$, 



(10) 



where /(f) = Y^f, M<l>)x<t>- 

We now turn our attention to semi-conservative repli- 
cation on a single fitness peak landscape. This case is 
more complicated than for a conservative system, since 
viability genes often exist on both strands in nature. 
Hence, if there exists a sequence 0o with fitness a, it 
stands to reason that the sequence cf) should have fitness 
a as well, effectively creating a double fitness peak land- 
scape (this assumption is by no means fundamental to 
the work). However, noting that x^, = x^, for all times, 
both by definition and by conservation in Equation (10), 
this difficulty can be sidestepped. As long as n is not 
too small, the area around each fitness peak can be lo- 
cally treated as a single fitness peak landscape as the 
two peaks are distant in sequence space. Hence, ignor- 
ing back mutations, the two master sequences obey the 
equations 



dwQ 
~df 



= 2(1 - e/2) n aw - (a + f(t))w = 



2(l-e/2) n aw -(a + f(t))w 



(11) 



dw 
~dt 



± = 2(1 - e/2) n aw - (a + f(t))w = 



2(1 -e/2) n aw - (a + /(f)K, 



(12) 



where Wi represents the concentration of the ith Ham- 
ming class as before. Therefore, we can re-define the 
concentration of the master sequence to include both wq 
and w and use equation (11) for the sum of the two. 
While this is not strictly necessary and has no effect on 
the results, it does reduce the bookkeeping, and the char- 
acteristics of the individual peaks can be obtained by 
simply dividing by two. A similar procedure yields 

dw-i , , K „ i , e/2 , 
^=2 ( l-e/2r-^)n™o + 



2(l-e/2)"r 7Wl -(7 ? + /(f)) Wl , 



(13) 
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where we include sequences of Hamming distance one 
away from both master sequences. The definition = 
Wiex.p(J* f(s)ds) once again removes the non-linearity. 
The solutions for the first two Hamming classes are 



i/o(t) = y (0)e 2 ^ 1 - £ / 2 ) r 



2/i (t) = Vo (0)n 



(14) 



ae(l -e/2) r 



(5-l)(a-r ? )(2(l-e/2)«-l) 



<r(2(l-e/2)"-l)t _ «(2(l-e/2)"-l)t 



(15) 



III. HOST-PARASITE CO-EVOLUTION 

Historically, the main focus of research on the quasis- 
pecies model has related to static and equilibrium prop- 
erties of the gystemSiSiiiiliiSiiSiii. A number of recent 
works, however, have explored the dynamics of the sys- 
tem under various conditionsiSiiSiiS, which has allowed 
the study of the simple model of co-evolution described 
here. Following the work of Kamp and BornholdliSiiS, 
we envision a population of host and parasite organisms 
(which we shall refer to as the immune system and virus) , 
each described by a set of quasispecies equations. Ignor- 
ing the interspecies interaction, the immune and viral 
genomes, of length n; s and n v , respectively, evolve inde- 
pendently on a single fitness peak landscape, where the 
master sequences have fitness Oi S rji a and <j v r/ v . 
To model the deleterious effect of the immune system on 
the virus, the dominant immune genome imposes a large 
death rate S on the corresponding viral sequence. If this 
dominant immune genome matches the viral master se- 
quence, the viral fitness peak will move to an arbitrary 
sequence of the first Hamming class. The viral quasis- 
pecies then adapts to this new fitness peak on a timescale 
t v , the time required for the population of the new mas- 
ter sequence to overtake that of the old. At this point, 
the immune system fitness peak adjusts to match the new 
viral peak, and adapts on a similarly defined timescale 
Ti s . Thus, through the iteration of these steps, the viral 
fitness peak scours sequence space in an attempt to avoid 
the immune system, which follows on its heels. Applying 
recent results on dynamic fitness landscapes^, regions of 
stability can be defined for both the viral and immune 
quasispecies by determining a characteristic timescale for 
regrowth of a new master sequence. If the landscape 
moves slowly enough, the master sequence has time to re- 
generate to the master sequence concentrations reached 
before the peak shift and the species will survive for all 
time. If, however, the master sequence cannot regener- 
ate rapidly enough, a second peak shift will occur before 
the new master sequence reaches the concentration held 
by the old master sequence before the first shift. The 
third master sequence cannot reach the levels of the sec- 
ond, and this continues until, eventually, there is no dis- 



cernible master sequence in the population. For the con- 
servative case, this can be stated rigorously by compar- 
ing the growth of a single member of the first Hamming 
class described by equation (8) with e 7?T , the uninhibited 
growth of a random sequence far from the fitness peak 
(as mutations in and out of this sequence should cancel). 
Using equation (8) this ratio can be defined, for both the 
immune and viral quasispecies, as&£ 

_ w\{t) 



(S-l)(<J- V )q 



(16) 



where r is the lag time between peak shifts and the pa- 
rameters {q, a, • • • } represent the parameters for either 
species. The quasispecies survives only when k > 1. 

The last piece necessary to complete the co-evolution 
model, then, is the speed with which the landscape 
moves. By the definition of our model, r is the sum 
of the the time required for the regeneration of the virus, 
t v , plus the time required for the regeneration of the im- 
mune system, Ti S . Hence, we must solve for r = t,- is + t v 
where 



This can be solved to obtain 



(17) 



(18) 



(S - l){a v - r] v )q v 



(19) 



jlruT* ( eq " a " T ~ g^Kl - q ia )<Tis (20) 
(S - l)(cr is - r] is )q is 

which yields, with the reasonable approximations that 
gsil and a ^ i] (the latter of which is used throughout 
the paper), 



(21) 



(22) 



These equations can be applied to determine the optimal 
mutation rate for both the host and the parasite. The 
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host can minimize the region of viability for the parasite 
by evolving a mutation rate such that 



de ls 



0. 



yielding^ 



1 - ni s ei s In 



S- 1 



(23) 



(24) 



of the parameters of the immune model, as well as the 
properties of the virus. The solution to this equation 
is shown in Fig. 2 and compared to the experimentally 
verified mutation rate for human B-cell receptors. This 
is discussed at length in section IV. 

Optimizing the viral mutation rate requires solving for 



dn v 



= 



(25) 



This equation has the nice quality of being independent o»Afi 



(Qv"( a v ~ Vv) + 5){n v (q v - l)ql n " alT is + 8(q v + (q v - l)n v q™ v a v T is ) + r v (q v - - (q v - l)n v q 2 v nv a v T is )) 



+n v q" v (q v - l)(rfi - 8a v - rj v a v )ln( g " ) = 0, 



(26) 



the solution of which is shown in Fig. 3 for a chosen set 
of parameters. 

We now turn our attention to the central theme of this 
paper, the co-evolution of semi-conservative organisms. 
Applying the results of section II and following the pro- 
cedure outlined above, we find, for a scmi-conscrvatively 
replicating host, 

a is e is {l-e is /2) n ^- x 



1> (2n v (l-e v /2) n "-2r, v )r 



(S-l)(a is -r) is )(2(l-e is /2y 



(2cr i( ,(l-ei«/2) r *"-(T is -r) is )r 



(2i7 is (l-ei s /2) n is_2 Ws ) 7 



'is i ' v 



1 ( (l-e ie /2) n t'e is \ 
m l(2(l- £ „/2)", s -l)(S-l)J 

(2(1 -£ is /2)"» -1)(«7. 



is Vis 



(27) 



(28) 



(29) 



A conservatively replicating virus interacting with this 
host will still follow the behavior described by equations 
16 and 21, albeit with the proper, semi-conservative r, s 
defined above. For the case of a semi-conservative virus 
we obtain 

<r„e„(l -e^/2)""- 1 



(S -l)(a v - J7B )(2(l-e,/2)» 



(e 



(2cr l ,(l-e„/2)"« -<J V -7] V )-T 



In 



(1-^/2)" 



(2(l-e„/2)«»-l)(S-l) 



(2(1 - e„/2)»» - l)(a v - r] v ) + S 
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Semi-conservative 



(30) 



(31) 




FIG. 2: Optimal immune system mutation rate vs. rii 3 . 
The dashed lines represent experimental values for somatic 
hypermutation of B-cell complementary determining regions, 
adapted from—. 



We now proceed to find the optimal mutation rates for 
both organisms. Differentiating k v by ei S and setting the 
result to zero gives us a criterion for the optimal immune 
mutation rate, 
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(5-1) ( 2-(l-^) "'") 

-2 + e ls + n ls e is - 2 (l - " -2 + e is + n is e is ln( ^ — y — - J -) 



(1 - 2 (1 - ^f) n ") (e ia - 2) e ls {<j is - Vls ) 



= 0. 



(32) 



This equation has all of the nice properties of Equation 
(24), defining an optimal mutation rate for any genome 
length, independent of the parameters of the system. The 
solution to this equation is plotted in Fig. 2, along with 
the conservative solution and the experimental range for 
observed rates per base pair per generation of somatic hy- 
permutation in the complementary determining regions 
(CDR's) found in B-cell antigen receptors. 




FIG. 3: Optimal viral mutation rate vs. n v for a conser- 
vative and semi-conservative virus interacting with a semi- 
conservative immune system. rn a = 100, <Ti s = a v — 
100, n is =r) v = l,5 = 200, e ls = 0.001. 

To maximize the stability of the viral quasispecies we 
set dn v /de v = as before. After a fair bit of work, we 
obtain an unwieldy expression omitted here in the inter- 
est of space 20 . The expression simplifies immensely in 
the limit <5 — * oo , the limit of an ideally efficient immune 
system. In this limit, 

Thy C v 

2(1 -2(1- f)"-) 2 

""T- (1 2(~i =£f " 1 = °- (33) 

The ideally efficient immune system is not an unreason- 
able approximation, as immune systems are highly effi- 
cient in destroying invaders once a suitable antibody is 
produced. The full expression as well as the above lim- 
iting form are dependent on both the parameters of the 



model and the properties of the immune system as in the 
conservative case. The solution of the full expression for 
a particular set of parameters is shown in Fig. 3. 



IV. RESULTS AND DISCUSSION 

Given the fundamental differences between semi- 
conservative and conservative modes of replication, the 
most striking aspect of Fig. 2 and 3 is the similar- 
ity between the conservative and semi-conservative op- 
timal mutation rates at high n , particularly for the viral 
species. This is most easily understood by noting that, 
as (1 — e/2)' 1 — > 1 for any semi-conservatively replicating 




FIG. 4: (1 — e/2) n for the optimal mutation rate of a semi- 
conservative immune system and virus. This parameter can 
be used as a measure of the "conservativeness" of a semi- 
conservative system. a is = a v — 100,rji S = rj v = 1,5 = 
200, e is = 0.001. 

organism, the probability that a mutation will be found 
in the original strands after replication vanishes. Hence, 
in this limit, semi-conservative and conservative replica- 
tion are expected to mimic each other. This parameter is 
shown in Fig. 4 for the optimal viral and immune muta- 
tion rates. Clearly, with the exception of small immune 
genomes, the conservative system can be used as a good 
approximation for semi-conservative replication. It is im- 
portant to note, however, that this knowledge could not 
have been extracted from the data for the conservative 
system. A large value for (1 — e/2)™ in the conserva- 
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tive system is a necessary but not sufficient criterion to 
justify the use of a conservative model, and the full semi- 
conservative calculation is required. 




10- 6 1CT 4 g i<r 2 10° 



FIG. O'. Kv VS. Ev for a conservative and semi-conservative 
virus interacting with a semi-conservative immune system. 
riis = n v — 100, ais = <J V = 100,?7j S = rj v = 1,5 = 200, tis = 
0.001. 

Equation (24) remains dependent on the parameters 
of the model, but general trends are obvious when bio- 
logically reasonable parameters are employed. While the 
extremal behavior of Equation (27) and (30) differs little 
from Equation (16) for genome lengths that are not too 
small, the behavior away from the maxima differs greatly. 
Fig. 5 displays n v vs. e v for a given set of parameters 
for both the conservative and semi-conservative models. 
It is immediately clear that, while the two models co- 
incide at small e (with a slightly higher peak height for 
either species for some parameters), their behavior dif- 
fers greatly otherwise, with the semi-conservative model 
displaying a more drastic dropoff in viability as e in- 
creases, true for all biologically reasonable parameters 
studied. The parameters shown in Fig. 5 were chosen 
as a representative, rather than extreme, example of this 
behavior. The importance of this result is best under- 
stood in light of the evolutionary pressures one would 
expect a viral population to encounter. The indepen- 
dence of Equation (32) from the properties of the viral 
system suggests that there exists an optimal mutation 
rate for an immune receptor independent of the qualities 
of the parasite against which it is defending. Thus, it 
is reasonable to expect (within the limitations imposed 
by additional evolutionary pressures, such as the need 
to distinguish between self and foreign antigens) an im- 
mune receptor to evolve this mutation rate nearly exactly. 
However, in the viral case, the optimal mutation rate de- 
pends strongly on the nature of the immune system it 
is attacking. Thus, the virus must evolve the mutation 
rate that maximizes its overall viability against the range 



of immune systems it is likely to infect, including both 
inter- and intra-species viability. The mutation rate that 
optimizes defense against one host may be a poor choice 
for another, and the virus must find the mutation rate 
that affords the best protection against all hosts, even if 
this is not the best mutation rate for evading any partic- 
ular immune system. Such a compromise clearly involves 
the behavior of k v over a wide range of e, rather than just 
at the maximum. One would therefore expect the more 
drastic dropoff at higher e to force the semi-conservative 
virus to develop a lower mutation rate so as to increase 
its viability against immune systems that lower the e„ 
with the maximal value of dn v /de v . Quantifying this 
statement requires an intelligent estimate of the distribu- 
tion of immune properties, a subject of future research. 
Qualitatively, this agrees well with the experimentally 
verified fact that semi-conservative viruses display sig- 
nificantly lower mutation rates than their conservative 
counterparts2ii22. 



V. CONCLUSIONS 

In this paper, we have extended Kamp and Born- 
holdt's model of co-evolution to incorporate the semi- 
conservative nature of DNA replication for both species. 
A parameter-independent expression was derived for the 
optimal mutation rate of an immune receptor, which 
agrees well with experimental data. Convergence of the 
conservative and semi-conservative results was demon- 
strated for realistic genome sizes, justifying the use of a 
conservative model in this case. 

Optimizing the stability of the immune species yielded 
a maximum that coincides with the conservative model 
for realistic genome sizes. A similar correspondence ex- 
ists for the virus, albeit with a dependence on the pa- 
rameters of the model. Away from the maximum, the 
conservative and semi-conservative models display differ- 
ent behaviors that provides a possible explanation for the 
high mutation rates found in conservative viruses. It is 
always dangerous to extrapolate from a simplified model 
of this kind to the complex systems found in nature. A 
true virus and immune system must contend with innu- 
merable evolutionary pressures, biological, chemical and 
otherwise, such as the requirement that T-cells recognize 
and do not bind host proteins. The work represented in 
this paper describes a generalized model which we feel 
captures the robust qualitative features of host-parasite 
coevolution, providing insight into the complex workings 
of nature. 
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